Failure avalanches in fiber bundles for discrete load increase 
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The statistics of burst avalanche sizes n during failure processes in a fiber bundle follows a power 
law, D{n) ~ for large avalanches. The exponent £ depends upon how the avalanches are 

provoked. While it is known that when the load on the bundle is increased in a continuous manner, 
the exponent takes the value £ = 5/2, we show that when the external load is increased in discrete 
and not too small steps, the exponent value £ = 3 is relevant. Our analytic treatment applies to 
bundles with a general probability distribution of the breakdown thresholds for the individual fibers. 
The pre-asymptotic size distribution of avalanches is also considered. 

PACS numbers: 62.20.Mk 



I. INTRODUCTION 

Fiber bundles with statistically distributed thresholds 
for breakdown of individual fibers provide simple and in- 
teresting models of failure processes in materials under 
stress. These models are much studied since they can 
be analyzed to an extent that is not possible for more 
complex materials (for reviews see [ITS! B 13). We 
study here equal-load-sharing models, in which the load 
previously carried by a failed fiber is shared equally by 
all the remaining intact fibers in the bundle 0, 0, \0, 3- 
The break-down thresholds Xi for the individual fibers 
are assumed to be independent random variables with 
the same cumulative distribution function P{x) and a 
corresponding probability density p(x): 

r 

Prob(iEi < x) = P(x) = I p(y) dy. (1) 
Jo 

We consider a bundle of N fibers, clamped at both 
ends. At a force x per surviving fiber the total force F(x) 
on the bundle is x times the number of intact fibers. The 
average, or macroscopic, force is given by the expectation 
value of this, 

(F) = N x [1 - P(x)]. (2) 

One may consider x to represent the elongation of the 
bundle, with the elasticity constant set equal to unity. 
In the generic case (F) will have a single maximum F C) a 
critical load corresponding to the maximum load the bun- 
dle can sustain before complete breakdown of the whole 
system. 

When a fiber ruptures somewhere, the stress on the 
intact fibers increases. This may in turn trigger further 
fiber failures, which can produce avalanches that either 
lead to a stable situation or to breakdown of the whole 
bundle. One may study the burst distribution D(n), de- 
fined as the expected number of bursts of size n when the 
bundle is stretched until complete breakdown. When the 
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load on the bundle is increased continuously from zero, 
the generic result is a power law 

Urn iV" 1 D(n) oc n~ s , (3) 

for large n, with £ = 5/2 (TSUI- 

However, experiments may be performed in a different 
manner. In Sec. 2 we show for the uniform probability 
distribution of thresholds that for experiments in which 
the load is increased in finite steps of size 5 rather than 
by infinitesimal amounts, the power-law exponent is seen 
to increase to the value £ = 3, This have been noticed 
in a special case previously: An argument by Pradhan et 
al. [12] for uniform threshold distribution suggested this 
exponent value. 

In Sec. 3 we show that the same asymptotic power law, 

D(n) oc 7tT 3 (4) 

holds for a triangular threshold distribution and for a 
Weibull distribution. The result |(4]) is demonstrated by 
simulations, and we provide analytic derivations to back 
up the results. Our conclusions are not limited to the 
distributions used in the simulations; we give an analytic 
derivation valid for a general threshold distribution. 

The avalanche distribution does not follow the asymp- 
totic power law |(4]) for small bursts. However, our ana- 
lytic formulas do also cover the distribution for bursts of 
smaller sizes n, down to a minimum size. 

In the concluding remarks we discuss briefly the results 
and their range of applicability. 



II. UNIFORM THRESHOLD DISTRIBUTION 

For the uniform distribution of thresholds, 

D , , / x for < x < 1 , , 

P ^ = \0 forrOl, ^ 

here in dimensionless units, the load curve ^ is 
parabolic, 

(F) = N x (1 - x), (6) 



2 



so that the expected critical load equals 

F c = N/4. (7) 

Fig. 1 shows the size distribution of the bursts obtained 
by simulation on fiber bundles with the uniform threshold 
distribution . In one procedure the bursts are recorded 
under a continuous load increase, in the other procedure 
the load is increased by discrete amounts 6. Clearly the 
exponent describing the large-size power laws are differ- 
ent, close to £ = 2.5 and £ = 3, respectively. 



the importance of force fluctuations [10]. An effective 
reduction of the fluctuations requires that the size 5 of 
the load increase is large enough so that a considerable 
number of fibers break in each step, i.e. 6 ^> F c /N. 

With a sufficiently large 6 we may use the macroscopic 
load equation to determine the number of fibers bro- 
ken in each step. The load values are mS, with m taking 
the values to = 0, 1, 2, ... , N/46 for the uniform threshold 
distribution. By ((6]) the threshold value corresponding to 
the load m6 is 
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FIG. 1. Avalanche size distribution for the uniform threshold 
distribution J5j. The distribution marked (a) is a record of 
all bursts when the load is increased continuously, while the 
bursts recorded in (b) results from increasing the load in steps 
of S = 10 and 8 = 50 (upper curve). The dotted line in (a) 
represents a n -5 / 2 behavior, and the dotted lines in (b) show 
the theoretical asymptotics l|12|l for S = 10 and <5 = 50. The 
figure is based on 10000 samples with N = 10 6 fibers in the bundle. 

The basic reason for the difference in the power laws 
is that increasing the external load in steps reduces the 
fluctuations in the force. The derivation of the asymp- 
totic size distribution D(n) oc n~ 5 / 2 of avalanches, corre- 
sponding to load increases by infinitesimal steps, shows 



4m6/NJ . 



(8) 



The expected number of fibers broken when the load is 
increased from m6 to (m + 1)6 is close to 



Ndx m /dm = S/y/l - 4m6/N. 



(9) 



Here the minimum number of n is 6, obtained in the first 
load increase. The integral over all m from to N/46 
yields a total number N/2 of broken fibers, as expected, 
since the remaining one-half of the fibers burst in one 
final avalanche. 

The number of avalanches of size between n and n+dn, 
D(n) dn, is given by the corresponding interval of the 
counting variable to: 



D(n) dn = dm. 



(10) 



Since 



dn _ 26 2 
dm~ N 



(1 - 4m6/N)~ 3 / 2 = — c n 3 
N6 



(11) 



we obtain the following distribution of avalanche sizes: 



D( n) = ^ = ±N6 n" 3 , 
dn z 



(n>6). (12) 



For consistency, one may estimate the total number of 
bursts by integrating D(n) from n — 6 to oo, with the 
result N/46, as expected. 

Fig. 1 shows that the theoretical power law (fT2"ll fits 
the simulation results perfectly for sufficiently large n. 
The simulation records also a few bursts of magnitude 
less than i5 because there is a nonzero probability to have 
bundles with considerably fewer fibers than the average 
in a threshold interval. However, these events will be 
of no importance for the asymptotic power law in the 
size distribution. Moreover, for a more realistic threshold 
distribution the situation is different (see below). 
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III. GENERAL THRESHOLD DISTRIBUTIONS 
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FIG. 2 Avalanche size distribution for (a) the triangular 
distribution lfl3ll and (b) for the Weibull distribution iflill . In both 
cases the load is increased in steps of S = 20. The figure is based 
on 10000 samples of bundles with N = 10 6 fibers. 

In order to see whether the asymptotic law (@| is gen- 
eral, we have performed simulations for two other thresh- 
old distributions, the triangular threshold distribution 



p(x) 



2x for < x < 1 
for x > 1, 



and the Weibull distribution with index 5, 

P(x) = 1 — e~ x , for all positive x. 



(13) 



(14) 



In Fig. 2 the exponent £ = 3 is clearly present in 
the triangular distribution case, and less clearly in the 



Weibull distribution case. Therefore we give now an an- 
alytic derivation for a general threshold distribution, and 
apply the result to the two cases for which the simula- 
tion results are shown in Fig. 2. The intention is to see 
if the simulation results can be fitted, not merely with a 
power law behavior for large avalanche sizes, but over an 
extended range of n. 

For a general threshold distribution P(x) a load inter- 
val 8 and a threshold interval are connected via the load 
equation ([2]). Since d(F)/dx = N[l — P(x) — xp(x)], an 
increase 8 in the load corresponds to an interval 



(15) 



N[l - P(x) - xp{x)] 



of fiber thresholds. The expected number of fibers broken 
by this load increase is therefore 



n = N p(x) Ax 



p(x) 



1 — P(x) — xp(x) 



8. 



(16) 



Note that this number diverges at the critical point, i.e. 
at the maximum of the load curve, as expected. 

We want to determine the number D(n) An of bursts 
with magnitudes in the interval (n, n + An). The number 
D{n) An we seek is the number of load steps correspond- 
ing to the interval An. Since each load step corresponds 
to the interval (fT5|) of x, the number of load steps corre- 
sponding to the interval An equals 



D(n) An = An 
Using (fl6|) we have 



dx N[l - P(x) - xp(x)) 
dn 8 



dn _ [1- P{x)] p'{x) + 2p{xf 
dx [1 — P(x) — xp(x)] 2 

Thus we obtain 

[1 - P{x) - xp(x)} 3 



(17) 



(18) 



D(n) 



N 

[1 - P{x)]p'{x) + 2p(x) 2 P 

p{xf N8 
[1 - P(x)]p'(x) + 2p{x) 2 n^"' 



(19) 



using (fT6|) . Here x = x(n), determined by Eq. (fT6|) . For 
the uniform threshold distribution expression (fl9|) coin- 
cides with (jl2| . 

Near criticality the first fraction in the expression (fT9|) 
for D(n) becomes a constant. The asymptotic behavior 
for large n is therefore 



D(n) 

with a nonzero constant 

C = N8 



Cn 



2p(x c ) + x c p'(x c ) 



(20) 



(21) 



We have used that at criticality 1 — P(x c ) = x c p(x c ). 
Thus the asymptotic power law (j4]) is universal. 
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For the triangular threshold distribution considered in 
Fig. 2a equations JH]) and CGI yield n = 2x6/(1 - 3x 2 ) 
and D = 4N6n~ 3 x 3 / (1 + 3a: 2 ). Elimination of x gives 



n 3 66 /n + 26 3 /n s + (3 + 25 2 /n 2 ) 1 /3 + S 2 /n 2 

(22) 

with asymptotic behavior D(n) ~ (2iV5/v27) n 3 , in 
agreement with (J20J) - 

For the Weibull distribution considered in Fig. 2b we 
obtain 

and n=T _ (23) 

This burst distribution must be given on parameter form, 
the elimination of x cannot be done explicitly. The criti- 
cal point is at x — and the asymptotics is given by 
(|2f)f. with C = N6(625e)~ 1 / 5 . 

These theoretical results are also exhibited in Fig. 2. 
In both cases the agreement with the theoretical results 
is very satisfactory. The n~ 3 power law is seen, but the 
asymptotics sets in only for very large avalanches, espe- 
cially for the Weibull case. For small-sized avalanches it 
is interesting to note the difference between the uniform 
distribution and the more realistic Weibull distribution. 
In the latter case the theoretical result (f!9l) is reasonably 
accurate also for n < 6. The reason is that for the Weibull 
(and the triangular) threshold distribution there are few 
very weak fibers. Thus a load increase 6 corresponds to 
a threshold interval that may, in the average description, 
contain just a few failing fibers. For the uniform dis- 
tribution on the other hand, a corresponding threshold 
interval with a number of fibers much less than 6 can only 
be caused by fluctuations with very low probabilities. 



IV. CONCLUDING REMARKS 

If we let the load increase 5 shrink to zero, we must 
recover the asymptotic D(n) oc n~ 5 / 2 power law valid for 
continuous load increase. Thus, as function of 5, there 
must be a crossover from one behavior to the other. It is 
to be expected that for 6 <C 1 the D(ri) oc ra~ 5 / 2 asymp- 
totics is seen, and when 5 > 1 the D(n) oc n~ 3 asymp- 
totics is seen. 

We would like to emphasize that the step-wise load 
increase is a reasonable as well as more practical load- 
ing method from the experimental point of view. While 
performing fracture-failure experiments by applying ex- 
ternal load, one cannot ensure that a single fiber (the 
weakest one) fails, whereas increasing the external load 
by equal steps is an easier and more realistic procedure. 

In conclusion we have shown that the magnitude dis- 
tribution D(n) of avalanches generated by simulation for 
fiber bundles with step- wise increase in the external load, 
is asymptotically a power law with an exponent essen- 
tially equal to —3. We have analytically derived this 
asymptotic power law D(n) = CrT 3 for a general proba- 
bility distribution p(x) of the individual fiber thresholds, 
as well as the pre-asymptotic behavior of the avalanche 
distribution D(n). 
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